Key mRNAs and lncRNAs of pituitary that affect the reproduction of FecB + + small tail han sheep.

BACKGROUND: The pituitary directly regulates the reproductive process through follicle-stimulating hormone (FSH) and luteinizing hormone (LH). Transcriptomic research on the pituitaries of ewes with different FecB (fecundity Booroola) genotypes has shown that some key genes and lncRNAs play an important role in pituitary function and sheep fecundity. Our previous study found that ewes with FecB + + genotypes (without FecB mutation) still had individuals with more than one offspring per birth. It is hoped to analyze this phenomenon from the perspective of the pituitary transcriptome. RESULTS: The 12 Small Tail Han Sheep were equally divided into polytocous sheep in the follicular phase (PF), polytocous sheep in the luteal phase (PL), monotocous sheep in the follicular phase (MF), and monotocous sheep in the luteal phase (ML). Pituitary tissues were collected after estrus synchronous treatment for transcriptomic analysis. A total of 384 differentially expressed genes (DEGs) (182 in PF vs. MF and 202 in PL vs. ML) and 844 differentially expressed lncRNAs (DELs) (427 in PF vs. MF and 417 in PL vs. ML) were obtained from the polytocous-monotocous comparison groups in the two phases. Functional enrichment analysis showed that the DEGs in the two phases were enriched in signaling pathways known to play an important role in sheep fecundity, such as calcium ion binding and cAMP signaling pathways. A total of 1322 target relationship pairs (551 pairs in PF vs. MF and 771 pairs in PL vs. ML) were obtained for the target genes prediction of DELs, of which 29 DEL-DEG target relationship pairs (nine pairs in PF vs. MF and twenty pairs in PL vs. ML). In addition, the competing endogenous RNA (ceRNA) networks were constructed to explore the regulatory relationships of DEGs, and some important regulatory relationship pairs were obtained. CONCLUSION: According to the analysis results, we hypothesized that the pituitary first receives steroid hormone signals from the ovary and uterus and that VAV3 (Vav Guanine Nucleotide Exchange Factor 3), GABRG1 (Gamma-Aminobutyric Acid A Receptor, Gamma 1), and FNDC1 (Fibronectin Type III Domain Containing 1) played an important role in this process. Subsequently, the reproductive process was regulated by gonadotropins, and IGFBP1 (Insulin-like Growth Factor Binding Protein 1) was directly involved in this process, ultimately affecting litter size. In addition, TGIF1 (Transforming Growth Factor-Beta-Induced Factor 1) and TMEFF2 (Transmembrane Protein With EGF Like And Two Follistatin Like Domains 2) compensated for the effect of the FecB mutation and function by acting on TGF-ß/SMAD signaling pathway, an important pathway for sheep reproduction. These results provided a reference for understanding the mechanism of multiple births in Small Tail Han Sheep without FecB mutation.

Comparative Transcriptomics Identify Key Pituitary Circular RNAs That Participate in Sheep (Ovis aries) Reproduction.

CircRNAs have been found to play key roles in many biological processes and have diverse biological functions. There have been studies on circRNAs in sheep pituitary, and some important circRNAs have been found. But there are still few studies on circRNAs in sheep pituitary with different fecundity. In this study, we obtained the circRNAs expression profiles in the pituitary of FecB ++ genotype Small Tail Han sheep with different fecundity and estrous phases. A total of 34,878 circRNAs were identified in 12 pituitary samples, 300 differentially expressed circRNAs (DE circRNAs) (down: 104; up: 196) were identified in polytocous sheep in the follicular phase (PF) and monotocous sheep in the follicular phase (MF) (PF vs. MF), and 347 DE circRNAs (down: 162; up: 185) were identified in polytocous sheep in the luteal phase (PL) and monotocous sheep in the luteal phase (ML) (PL vs. ML). Cortisol synthesis and secretion pathway (follicular phase) and estrogen signaling pathway (luteal phase) were obtained by functional enrichment analysis of circRNAs source genes. Competing endogenous RNA (ceRNA) network analysis of key DE circRNAs revealed that oar-circ-0022776 (source gene ITPR2, follicular phase) targeted oar-miR-432, oar-circ-0009003 (source gene ITPR1, luteal phase) and oar-circ-0003113 (source gene PLCB1, luteal phase) targeted oar-miR-370-3p. We also explored the coding ability of DE circRNAs. In conclusion, our study shows that changes in the pituitary circRNAs may be related to the response of the pituitary to steroid hormones and regulate the reproductive process of sheep by affecting the pituitary function. Results of this study provide some new information for understanding the functions of circRNAs and the fecundity of FecB ++ genotype sheep.

Bovine colostrum promoted ileal health in newborn lambs at 24 h after birth: insight from intestinal morphology and innate immunity.

The contribution of colostrum to passive immunity transfer and intestinal protection in newborn ruminants is well known; however, it is currently unclear how colostrum intake affects intestinal innate immunity. We investigated the effects of bovine colostrum intake on ileal morphology, expression of genes involved in intestinal innate immunity, and serum concentrations of inflammatory cytokines in newborn lambs. Twenty-seven newborn male Hu lambs were used, of which 18 were bottle-fed either bovine colostrum (C24h; n = 9) or bovine mature milk (M24h; n = 9) within the first 2 h after birth at an intake of approximately 8% of BW; the remaining nine lambs did not receive any feeding (N24h). Blood and ileal tissue samples were collected after the lambs were slaughtered at 24 h after birth. Ileal villus height and villus height-to-crypt depth ratio were significantly higher in C24h than those in N24h and M24h lambs (P < 0.01). Messenger RNA (mRNA) abundance of toll-like receptor (TLR)-2, TLR3, TLR4, TLR6, TLR7, TLR8 and tumour necrosis factor alpha in the ileum was lower in C24h than that in N24h lambs (P < 0.05). Moreover, C24h lambs had a lower TLR3 mRNA abundance (P < 0.01) and a trend of lower TLR6 (P = 0.06) and interleukin 1 beta (P = 0.08) expression compared with those in M24h lambs. We also observed strong positive correlations of tumour necrosis factor alpha expression with that of TLR2 (r = 0.71; P < 0.001), TLR4 (r = 0.88; P < 0.001) and TLR8 (r = 0.83; P < 0.001). Interestingly, the expression of barrier-related molecules, including mucin-13, lysozyme, claudin (CLDN)-1, CLDN2, CLDN4, CLDN7, CLDN12, occludin, zonula occluden-1 and junctional adhesion molecule-1, was consistently lower in C24h lambs than that in N24h and M24h lambs (P < 0.05). These results indicated that the beneficial roles of colostrum intake on intestinal protection in newborn lambs were associated with low TLR expression, which was reflected by improved intestinal development and reduced inflammatory response. Further studies using fluorescence in situ hybridisation and immunohistochemical methods are needed to further explore the mechanisms underlying the lower expression of intestinal barrier-related molecules due to colostrum feeding.

Expression, structure and function analysis of the sperm-oocyte fusion genes Juno and Izumo1 in sheep (Ovis aries).

BACKGROUND: JUNO and IZUMO1 are the first receptor-ligand protein pairs discovered to be essential for sperm-oocyte fusion; their interaction is indispensable for fertilization. METHODS: PCR was used to clone the full-length DNA sequence of the Juno gene in sheep. The single nucleotide polymorphism (SNP) loci of Juno were genotyped by Sequenom MassARRAY®. PCR combined with rapid amplification of cDNA Ends were used to clone the full-length cDNA sequence of Juno and Izumo1. Reverse transcriptase-PCR (RT-PCR) and real time-quantitative-PCR (RT-qPCR) were used to analyze the genes' expression in tissues of sheep, and single cell RNA-seq was used to analyze the genes' expression in oocytes, granulosa cells and follicular theca of polytocous and monotocous Small Tail Han ewes. Bioinformatics was used to analyze advanced structure and phylogeny of JUNO and IZUMO1 proteins. RESULTS: The full-length DNA sequence of the Juno gene in sheep was cloned and nine SNPs were screened. We found a significant association between the g.848253 C > A locus of Juno and litter size of Small Tail Han sheep (P < 0.05). The full-length cDNA sequence of Juno and Izumo1 genes from Small Tail Han sheep were obtained. We found a new segment of the Izumo1 CDS consisting of 35 bp, and we confirmed the Izumo1 gene has 9 exons, not 8. RT-qPCR showed that Juno and Izumo1 genes were highly expressed in ovarian and testicular tissues, respectively (P < 0.01). Single cell RNA-seq showed Juno was specifically expressed in oocytes, but not in granulosa cells or follicular theca, while Izumo1 displayed little to no expression in all three cell types. There was no difference in expression of the Juno gene in oocyte and ovarian tissue in sheep with different litter sizes, indicating expression of Juno is not related to litter size traits. Bioinformatic analysis revealed the g.848253 C > A locus of Juno results in a nonconservative missense point mutation leading to a change from Phe to Leu at position 219 in the amino acid sequence. CONCLUSIONS: For the first time, this study systematically analyzed the expression, structure and function of Juno and Izumo1 genes and their encoded proteins in Small Tail Han sheep, providing the basis for future studies of the regulatory mechanisms of Juno and Izumo1 genes.

Polymorphism, expression and structure analysis of key genes in the ovarian steroidogenesis pathway in sheep (Ovis aries).

BACKGROUND: Litter size is an important factor that significantly affects the development of the sheep industry. Our previous TMT proteomics analysis found that three key proteins in the ovarian steroidogenesis pathway, STAR, HSD3B1, and CYP11A1, may affect the litter size trait of Small Tail Han sheep. OBJECTIVE: The purpose of this study was to better understand the relationship between polymorphisms of these three genes and litter size. MATERIAL AND METHOD: Sequenom MassARRAY detected genetic variance of the three genes in 768 sheep. Real-time qPCR of the three genes was used to compare their expression in monotocous and polytocous sheep in relevant tissues. Finally, bioinformatics analysis predicted the protein sequences of the different SNP variants. RESULT: Association analysis showed that there was a significant difference in litter size among the genotypes at two loci of the CYP11A1 gene (p < 0.05), but no significant difference was observed in litter size among all genotypes at all loci of the STAR and HSD3B1 genes (p > 0.05). However, STAR expression was significantly different in polytocous and monotocous sheep in the pituitary (p < 0.01). Tissue-specific expression in the ovary was observed for HSD3B1 (p < 0.05), but its expression was not different between polytocous and monotocous sheep. Bioinformatics analysis showed that the g.33217408C > T mutation of CYP11A1 resulted in major changes to the secondary and tertiary structures. In contrast, gene polymorphisms in STAR and HSD3B1 had minimal impacts on their protein structures. DISCUSSION: This may explain why the CYP11A1 variant impacted litter size while the others did not. The single nucleotide polymorphism of the CYP11A1 gene would serve as a good molecular marker when breeding to increase litter size in sheep. Our study provides a basis for further revealing the function of the ovarian steroidogenesis pathway in sheep reproduction and sheep breeding.

Identification and Characterization of Hypothalamic Alternative Splicing Events and Variants in Ovine Fecundity-Related Genes.

Previous studies revealed that alternative splicing (AS) events and gene variants played key roles in reproduction; however, their location and distribution in hypothalamic fecundity-related genes in sheep without the FecB mutation remain largely unknown. Therefore, in this study, we described the hypothalamic AS events and variants in differentially expressed genes (DEGs) in Small-tailed Han sheep without the FecB mutation at polytocous sheep in the follicular phase vs. monotocous sheep in the follicular phase (PF vs. MF) and polytocous sheep in the luteal phase vs. monotocous sheep in the luteal phase (PL vs. ML) via an RNA-seq study for the first time. We found 39 DEGs with AS events (AS DEGs) in PF vs. MF, while 42 AS DEGs were identified in PL vs. ML. No DEGs with single nucleotide polymorphisms (SNPs) were observed in PF vs. MF, but five were identified in PL vs. ML. We also performed a correlation analysis of transcriptomics and proteomics, and the results suggested several key DEGs/differentially expressed proteins (DEPs), such as LGALS3 in PF vs. MF and aspartoacylase (ASPA) and transthyretin (TTR) in PL vs. ML, could be candidate genes influencing ovine litter size. In addition, further analyses suggested that AS events, SNPs and miRNA-binding sites existed in key DEGs/DEPs, such as ASPA and TTR. All in all, this study provides a new insight into ovine and even other mammalian reproduction.

TGIF1 and SF1 polymorphisms are associated with litter size in Small Tail Han sheep.

TGF-ß induced factor homeobox 1 (TGIF1) and splicing factor 1 (SF1) are important for mammalian reproduction; however, the effects of these genes on litter size in sheep remain unexplored. In this study, we genotyped 768 ewes from seven sheep breeds at two loci: g.37871539C>T, a synonymous mutation of TGIF1; and g.42314637T>C, a 3'UTR variant of SF1. Our analysis of polymorphism revealed only two genotypes at locus g.37871539C>T in TGIF1, with most sheep populations being moderately polymorphic (0.25 < PIC < 0.5) at this site. In contrast, most breeds exhibited low polymorphism (PIC ≤0.25) at the SF1 locus g.42314637T>C. The association analysis revealed that a synonymous mutation at g.37871539C>T in TGIF1 was highly associated with litter size in Small Tail Han sheep, in which it causes a significant decrease in litter size. Conversely, while the SF1 3'UTR variant g.42314637T>C was also highly associated with litter size in sheep, it causes a significant increase in the number of litter size. Combined, these data provide valuable information regarding candidate genetic markers for sheep breeding programs.

Polymorphisms in the ASMT and ADAMTS1 gene may increase litter size in goats.

Prolificacy of most local goat breeds in China is low. Jining Grey goat is one of the most prolific goat breeds in China, it is an important goat breed for the rural economy. ASMT (acetylserotonin O-methyltransferase) and ADAMTS1 (ADAM metallopeptidase with thrombospondin type 1 motif) are essential for animal reproduction. Single nucleotide polymorphisms (SNPs) of ASMT and ADAMTS1 genes in the highly prolific breed (Jining Grey goats), medium prolific breed (Boer goats and Guizhou White goats) and low prolific breeds (Angora goats, Liaoning Cashmere goats and Inner Mongolia Cashmere goats) were detected by polymerase chain reaction-restriction fragment length polymorphism and sequencing. Two SNPs (g.158122T>C, g.158700G>A) of ASMT gene and two SNPs (g.7979798A>G, g.7979477C>T) of ADAMTS1 gene were identified. For g.158122T>C of ASMT gene, further analysis revealed that genotype TC or CC had 0.66 (p < 0.05) or 0.75 (p < 0.05) kids more than those with genotype TT in Jining Grey goats. No significant difference (p > 0.05) was found in litter size between TC and CC genotypes. The SNP (g.158122T>C) caused a p.Tyr298His change and this SNP mutation resulted in changes in protein binding sites and macromolecule-binding sites. The improvement in reproductive performance may be due to changes in the structure of ASMT protein. For g.7979477C>T of ADAMTS1 gene, Jining Grey does with genotype CT or TT had 0.82 (p < 0.05) or 0.86 (p < 0.05) more kids than those with genotype CC. No significant difference (p > 0.05) was found in litter size between CT or TT genotypes. These results preliminarily indicated that C allele (g.158122T>C) of ASMT gene and T allele (g.7979477C>T) of ADAMTS1 gene are potential molecular markers which could improve litter size of Jining Grey goats and be used in goat breeding.

Mutations in NLRP5 and NLRP9 are Associated with Litter Size in Small Tail Han Sheep.

Previous studies showed that the NLR family pyrin domain-containing 5 (NLRP5) and NLRP9 genes are two important reproductive genes; however, their effects on sheep litter size are unknown. Therefore, in this study, we first genotyped seven sheep breeds via the MassARRAY® SNP system at the loci g.60495375A > G, g.60495363G > A, and g.60499690C > A in NLRP5, and g.59030623T > C and g.59043397A > C in NLRP9. Our results revealed that each locus in most sheep breeds contained three genotypes. Then, we conducted population genetic analysis of single nucleotide polymorphisms in NLRP5 and NLRP9, and we found that the polymorphism information content value in all sheep breeds ranged from 0 to 0.36, and most sheep breeds were under Hardy-Weinberg equilibrium (p > 0.05). Furthermore, association analysis in Small Tail Han sheep indicated that two loci, g.60495363G > A in NLRP5 and g.59030623T > C in NLRP9, were highly associated with litter size. The mutation in g.60495363G > A may decrease interactions of NLRP5 with proteins, such as GDF9, whereas the mutation in g.59030623T > C may enhance the combining capacity of NLRP9 with these proteins; consequently, these mutations may influence the ovulation rate and even litter size. The findings of our study provide valuable genetic markers that can be used to improve the breeding of sheep and even other mammals.

Low Expression of Sirtuin 1 in the Dairy Cows with Mild Fatty Liver Alters Hepatic Lipid Metabolism.

Dairy cows usually experience negative energy balance coupled with an increased incidence of fatty liver during the periparturient period. The purpose of this study was to investigate the effect of hepatic steatosis on the expression of the sirtuin 1 (SIRT1), along with the target mRNA and protein expressions and activities related to lipid metabolism in liver tissue. Control cows (n = 6, parity 3.0 ± 2.0, milk production 28 ± 7 kg/d) and mild fatty liver cows (n = 6, parity 2.3 ± 1.5, milk production 20 ± 6 kg/d) were retrospectively selected based on liver triglycerides (TG) content (% wet liver). Compared with the control group, fatty liver cows had greater concentrations of cholesterol and TG along with the typically vacuolated appearance and greater lipid droplets in the liver. Furthermore, fatty liver cows had greater mRNA and protein abundance related to hepatic lipid synthesis proteins sterol regulatory element binding proteins (SREBP-1c), long-chain acyl-CoA synthetase (ACSL), acyl-CoA carbrolase (ACC) and fatty acid synthase (FAS) and lipid transport proteins Liver fatty acid binding protein (L-FABP), apolipoprotein E (ApoE), low density lipoprotein receptor (LDLR) and microsomal TG transfer protein (MTTP) (p < 0.05). However, they had lower mRNA and protein abundance associated with fatty acid ß-oxidation proteins SIRT1, peroxisome proliferator-activated receptor co-activator-1 (PGC-1α), peroxisome proliferator-activated receptor-α (PPARα), retinoid X receptor (RXRα), acyl-CoA 1 (ACO), carnitine palmitoyltransferase 1 (CPT1), carnitine palmitoyltransferase 2 (CPT2) and long- and medium-chain 3-hydroxyacyl-CoA dehydrogenases (LCAD) (p < 0.05). Additionally, mRNA abundance and enzyme activity of enzymes copper/zinc superoxide dismutase (Cu/Zn SOD), catalase (CAT), glutathione peroxidase (GSH-Px) and manganese superoxide dismutase (Mn SOD) decreased and mRNA and protein abundance of p45 nuclear factor-erythroid 2 (p45 NF-E2)-related factor 1 (Nrf1), mitochondrial transcription factor A (TFAM) decreased (p < 0.05). Lower enzyme activities of SIRT1, PGC-1α, Cu/Zn SOD, CAT, GSH-Px, SREBP-1c and Mn SOD (p < 0.05) and concentration of reactive oxygen species (ROS) were observed in dairy cows with fatty liver. These results demonstrate that decreased SIRT1 associated with hepatic steatosis promotes hepatic fatty acid synthesis and inhibits fatty acid ß-oxidation. Hence, SIRT1 may represent a novel therapeutic target for the treatment of the fatty liver disease in dairy cows.

Proteomic analysis of sheep uterus reveals its role in prolificacy.

Small Tail Han sheep have attracted attention for their high fecundity and year-round estrus. However, the molecular mechanisms of this fecundity are unknown. Polymorphism of the FecB gene has been shown to be associated with the ovulation rate and litter size in sheep. In this study, we used tandem mass tag quantitative proteomic techniques to identify the differentially abundant proteins in polytocous and monotocous Small Tail Han sheep (FecB++) uterine tissues in the follicular and luteal phases. In total, 41 and 43 differentially abundant proteins were identified in the follicular and luteal phases, respectively. Correlation analysis between the transcriptome and proteome revealed a positive correlation at the two omics levels of prolificacy. GO and KEGG pathway analyses indicated that the mRNAs and proteins upregulated in the polytocous group relative to the monotocous group are involved in sphingolipid metabolism and amino acid metabolism, and may be important in maintaining uterine functions and increasing the embryo survival rate during the estrus cycle of polytocous sheep. In conclusion, our work provides a prospective understanding of the molecular mechanism underlying the high prolificacy of Small Tail Han sheep. SIGNIFICANCE: Fecundity critically affects the profitability of sheep production, but the genetic mechanism of high-prolificacy is still unclear in sheep. We identified potential signaling pathways and differentially abundant proteins associated with reproductive performance through a combination of sheep uterus tissues proteome and transcriptome analyses. These findings will facilitate a better revealing the mechanism and provide possible targets for molecular design breeding for the formation of polytocous traits in sheep.

Differential Expression of Circular RNAs in Polytocous and Monotocous Uterus during the Reproductive Cycle of Sheep.

CircRNA plays important roles in cell proliferation, differentiation, autophagy and apoptosis during development. However, there are few reports on circRNAs related to livestock reproduction. In this study, we identified circRNAs by deep sequencing and analyzed their expression in the uteri of polytocous and monotocous sheep (FecB++) during follicular and luteal phases. There were 147 and 364 circRNAs with differential expression in the follicular and luteal phases, respectively. GO and KEGG enrichment analysis was performed for the host genes of the circRNAs to predict the functions of differentially expressed circRNAs. These source genes were mainly involved in the estrogen signaling pathway, TGFß signaling pathway, GnRH signaling pathway, oxytocin signaling pathway, pentose phosphate pathway, and starch and sucrose metabolism related to reproduction and energy metabolism. CircRNA expression patterns were validated by RT-qPCR. Our findings provide a solid foundation for the identification and characterization of key important circRNAs involved in reproduction.

Comparative Transcriptomics Identify Key Hypothalamic Circular RNAs that Participate in Sheep (Ovis aries) Reproduction.

Circular RNA (circRNA), as an emerging class of noncoding RNA, has been found to play key roles in many biological processes. However, its expression profile in the hypothalamus, a powerful organ initiating the reproductive process, has not yet been explored. Therefore, we used RNA sequencing to explore the expression of circRNAs in the hypothalamus of sheep with the FecB ++ genotype. We totally identified 41,863 circRNAs from sheep hypothalamus, in which 333 (162 were upregulated, while 171 were downregulated) were differentially expressed in polytocous sheep in the follicular phase versus monotocous sheep in the follicular phase (PF vs. MF), moreover, 340 circRNAs (163 were upregulated, while 177 were downregulated) were differentially expressed in polytocous sheep in the luteal phase versus monotocous sheep in the luteal sheep (PL vs. ML). We also identified several key circRNAs including oar_circ_0018794, oar_circ_0008291, oar_circ_0015119, oar_circ_0012801, oar_circ_0010234, and oar_circ_0013788 through functional enrichment analysis and oar_circ_0012110 through a competing endogenous RNA network, most of which may participate in reproduction by influencing gonadotropin-releasing hormone (GnRH) activities or affecting key gene expression, indirectly or directly. Our study explored the overall expression profile of circRNAs in sheep hypothalamus, which potentially provides an alternative insight into the mechanism of sheep prolificacy without the effects of FecB mutation.

Identification and characterization of mRNAs and lncRNAs in the uterus of polytocous and monotocous Small Tail Han sheep (Ovis aries).

BACKGROUND: Long non-coding RNAs (lncRNAs) regulate endometrial secretion and uterine volume. However, there is little research on the role of lncRNAs in the uterus of Small Tail Han sheep (FecB++). Herein, RNA-seq was used to comparatively analyze gene expression profiles of uterine tissue between polytocous and monotocous sheep (FecB++) in follicular and luteal phases. METHODS: To identify lncRNA and mRNA expressed in the uterus, the expression of lncRNA and mRNA in the uterus of Small Tail Han sheep (FecB++) from the polytocous group (n = 6) and the monotocous group (n = 6) using RNA-sequencing and real-time polymerase chain reaction (RT-PCR). Identification of differentially expressed lncRNAs and mRNAs were performed between the two groups and two phases . Gene ontology (GO) and pathway enrichment analyses were performed to analyze the biological functions and pathways for the differentially expressed mRNAs. LncRNA-mRNA co-expression network was constructed to further analyses the function of related genes. RESULTS: In the follicular phase, 473 lncRNAs and 166 mRNAs were differentially expressed in polytocous and monotocous sheep; in the luteal phase, 967 lncRNAs and 505 mRNAs were differentially expressed in polytocous and monotocous sheep. GO and KEGG enrichment analysis showed that the differentially expressed lncRNAs and their target genes are mainly involved in ovarian steroidogenesis, retinol metabolism, the oxytocin signaling pathway, steroid hormone biosynthesis, and the Foxo signaling pathway. Key lncRNAs may regulate reproduction by regulating genes involved in these signaling pathways and biological processes. Specifically, UGT1A1, LHB, TGFB1, TAB1, and RHOA, which are targeted by MSTRG.134747, MSTRG.82376, MSTRG.134749, MSTRG.134751, and MSTRG.134746, may play key regulatory roles. These results offer insight into molecular mechanisms underlying sheep prolificacy.

Identification of Prolificacy-Related Differentially Expressed Proteins from Sheep (Ovis aries) Hypothalamus by Comparative Proteomics.

Reproduction, as a physiologically complex process, can significantly affect the development of the sheep industry. However, a lack of overall understanding to sheep fecundity has long blocked the progress in sheep breeding and husbandry. In the present study, the aim is to identify differentially expressed proteins (DEPs) from hypothalamus in sheep without FecB mutation in two comparison groups: polytocous (PF) versus monotocous (MF) sheep at follicular phase and polytocous (PL) versus monotocous (ML) sheep at luteal phase. Totally 5058 proteins are identified in sheep hypothalamus, where 22 in PF versus MF, and 39 proteins in PL versus ML are differentially expressed, respectively. A functional analysis is then conducted including Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis to reveal the potential roles of these DEPs. The proteins ENSOARP00000020097, ENSOARP00000006714, growth hormone (GH), histone deacetylase 4 (HDAC4), and 5'-3' exoribonuclease 2 (XRN2) in PF versus MF, and bcl-2-associated athanogene 4 (BAG4), insulin-like growth factor-1 receptor (IGF1R), hydroxysteroid 11-beta dehydrogenase 1 (HSD11B1), and transthyretin (TTR) in PL versus ML appear to modulate reproduction, presumably by influencing the activities of gonadotropin-releasing hormone (GnRH). This study provides an alternative method to identify DEPs associated with sheep prolificacy from the hypothalamus. The mass spectrometry data are available via ProteomeXchange with identifier PXD013822.

The genetic mechanism of high prolificacy in small tail han sheep by comparative proteomics of ovaries in the follicular and luteal stages.

To investigate the genetic mechanism of sheep prolificacy, protein profiling of ovaries in the follicular and luteal phases was conducted. The tandem mass tag technique was used to analyze the proteomes of ovaries from STH sheep that did not have the FecB mutation in the bone morphogenetic protein receptor 1B gene. Parallel Reaction Monitoring (PRM) was operated to validate the target differentially abundant proteins (DAPs). The result showed, a total of 34,037 peptides were found, and 5074 proteins were identified. The screened DAPs strictly related to energy metabolism, hormone synthesis, ovarian function were significantly enriched in oxidative phosphorylation(COX7A, ND5, and UQCR10), ovarian steroidogenesis(StAR and HSD3B), taurine and hypotaurine metabolism(CSAD), glycosaminoglycan biosynthesis-heparin sulfate/heparin(GLCE), necroptosis(H2AX, AIFM1, and FTH1), protein digestion and absorption(COL4A1 and COL4A5) and glycosaminoglycan degradation(HYAL2 and HEXB) pathways. These analyses indicated that the reproductive performance of sheep is regulated through different pathways. In consequence, these findings are an important resource that can be used in future studies of the genetic mechanism of high fecundity traits in sheep, and these DAPs can be further investigated as candidate markers to predict prolificacy of sheep. SIGNIFICANCE: Litter size is an important quantitative trait, but the genetic mechanism of high-prolificacy is still unclear in sheep. Our study identified potential signaling pathways and differentially abundant proteins related to reproductive performance. These findings will facilitate a better revealing the mechanism and provide possible targets for molecular design breeding for the formation of polytocous traits in sheep.

Comparative Transcriptomics Reveal Key Sheep (Ovis aries) Hypothalamus LncRNAs that Affect Reproduction.

The diverse functions of long noncoding RNAs (lncRNAs), which execute their functions mainly through modulating the activities of their target genes, have been have been widely studied for many years (including a number of studies involving lncRNAs in the ovary and uterus). Herein, for the first time, we detect lncRNAs in sheep hypothalami with FecB++ through RNA Sequencing (RNA-Seq) and identify a number of known and novel lncRNAs, with 622 and 809 found to be differentially expressed in polytocous sheep in the follicular phase (PF) vs. monotocous sheep in the follicular phase (MF) and polytocous sheep in the luteal phase (PL) vs. monotocous sheep in the luteal phase (ML), respectively. Then, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed based on the predicted target genes. The most highly enriched GO terms (at the molecular function level) included carbonyl reductase (NADPH), 15-hydroxyprostaglandin dehydrogenase (NADP+), and prostaglandin-E2 9-reductase activity in PF vs. MF, and phosphatidylinositol-3,5-bisphosphate binding in PL vs. ML was associated with sheep fecundity. Interestingly, the phenomena of valine, leucine, and isoleucine degradation in PL vs. ML, and valine, leucine, and isoleucine biosynthesis in PF vs. MF, were present. In addition, the interactome of lncRNA and its targets showed that MSTRG.26777 and its cis-targets ENSOARG00000013744, ENSOARG00000013700, and ENSOARG00000013777, and MSTRG.105228 and its target WNT7A may participate in the sheep reproductive process at the hypothalamus level. Significantly, MSTRG.95128 and its cis-target Forkhead box L1 (FOXG1) were shown to be upregulated in PF vs. MF but downregulated in PL vs. ML. All of these results may be attributed to discoveries of new candidate genes and pathways related to sheep reproduction, and they may provide new views for understanding sheep reproduction without the effects of the FecB mutation.

Integrated Hypothalamic Transcriptome Profiling Reveals the Reproductive Roles of mRNAs and miRNAs in Sheep.

Early studies have provided a wealth of information on the functions of microRNAs (miRNAs). However, less is known regarding their functions in the hypothalamus involved in sheep reproduction. To explore the potential roles of hypothalamic messenger RNAs (mRNAs) and miRNAs in sheep without FecB mutation, in total, 172 and 235 differentially expressed genes (DEGs) and 42 and 79 differentially expressed miRNAs (DE miRNAs) were identified in polytocous sheep in the follicular phase versus monotocous sheep in the follicular phase (PF vs. MF) and polytocous sheep in the luteal phase versus monotocous sheep in the luteal phase (PL vs. ML), respectively, using RNA sequencing. We also identified several key mRNAs (e.g., POMC, GNRH1, PRL, GH, TRH, and TTR) and mRNA-miRNAs pairs (e.g., TRH co-regulated by oar-miR-379-5p, oar-miR-30b, oar-miR-152, oar-miR-495-3p, oar-miR-143, oar-miR-106b, oar-miR-218a, oar-miR-148a, and PRL regulated by oar-miR-432) through functional enrichment analysis, and the identified mRNAs and miRNAs may function, conceivably, by influencing gonadotropin-releasing hormone (GnRH) activities and nerve cell survival associated with reproductive hormone release via direct and indirect ways. This study represents an integral analysis between mRNAs and miRNAs in sheep hypothalamus and provides a valuable resource for elucidating sheep prolificacy.

Expression Analysis of the Prolific Candidate Genes, BMPR1B, BMP15, and GDF9 in Small Tail Han Ewes with Three Fecundity (FecB Gene) Genotypes.

The expression characteristics of the prolific candidate genes, BMPR1B, BMP15, and GDF9, in the major visceral organs and hypothalamic⁻pituitary⁻gonadal (HPG) axis tissues of three FecB genotypes (FecB BB, FecB B+, and FecB ++) were explored in STH ewes using RT-PCR and qPCR. The results were as follows, BMPR1B was expressed in all FecB BB genotype (Han BB) tissues, and GDF9 was expressed in all selected tissues, but BMP15 was specifically expressed in the ovaries. Further study of ovarian expression indicated that there was no difference in BMPR1B expression between genotypes, but the FecB B+ genotype (Han B+) had greater expression of GDF9 and BMP15 than Han BB and FecB ++ genotype (Han ++) (p < 0.05, p < 0.01). BMP15 expression was lower in the ovaries of Han BB than in Han ++ sheep, but the reverse was shown for GDF9. The gene expression in non-ovarian tissues was also different between genotypes. Therefore, we consider that the three genes have an important function in ovine follicular development and maturation. This is the first systematic analysis of the tissue expression pattern of BMPR1B, BMP15, and GDF9 genes in STH sheep of the three FecB genotypes. These results contribute to the understanding of the molecular regulatory mechanism for ovine reproduction.

Effects of deoxynivalenol exposure on cerebral lipid peroxidation, neurotransmitter and calcium homeostasis of chicks in vivo.

During current research, the effects of deoxynivalenol (DON) exposure on cerebral lipid peroxidation, neurotransmitter secretion and calcium homeostasis in chicks were evaluated. One hundred and twenty Hailan chicks (male, 1-day-old) were randomly divided into four groups. Chicks in low, medium and high dose groups were fed with 0.27, 1.68 and 12.21 mg/kg-1 DON respectively by gavage according to feed intake. Chicks in control group were fed with physiological saline by gavage. The trials were conducted for 36 d. At the end of the trials, twenty chicks per group were sacrificed, and the cerebra were collected for measuring the brain indices. Compared with the control group, the activities of total superoxide dismutase (T-SOD) and glutathione peroxidase were significantly decreased in treatment groups (P < 0.05), the contents of malondialdehyde in high dose group were increased (P < 0.05), the catalase activities and nitric oxide contents in medium and high dose groups were decreased (P < 0.05), and the activities of T-AOC in high dose group were reduced (P < 0.05). Compared with the control group, the concentrations of norepinephrine and 5-hydroxytryptamine in high dose group were obviously increased (P < 0.05), while the concentrations of dopamine were decreased (P < 0.05). Meanwhile, the concentrations of calcium and calmodulin (CaM) in medium and high dose groups were lower than those of the control group (P < 0.05), and the gene relative expression of CaM mRNA in treatment groups were significantly reduced (P < 0.05), in a dose-dependent manner. These results suggested that DON exposure can affect the cerebral lipid peroxidation, neurotransmitters secretion and the balance of calcium homeostasis in chicks.